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First-order phase transitions typically exhibit a significant hysteresis resulting for instance in boiling 
retardation and supercooling. The hysteresis arises, because nucleation of the new phase is activated. The 
free- energy change is positive until the nucleus reaches a critical size beyond which further growth is 
downhill. In practice, the barrier is often circumvented by the presence of heterogeneous nucleation centres, 
e.g. at vessel walls or seed crystals. Recently, it has been proposed that the homogeneous melting of ice 
proceeds via separation of defect pairs with a substantially smaller barrier as compared to the mere 
aggregation of defects. Here we report the observation of an analogous mechanism catalysing a 
two-dimensional homogeneous phase transition. A similar process is believed to occur in spin systems. This 
suggests that separation of defect pairs is a common trigger for phase transitions. Partially circumventing 
the activation barrier it reduces the hysteresis and may promote fluctuations within a temperature range 
increasing with decreasing dimensionality. 

Phase transitions are well understood from a thermodynamic point of view, but the kinetics of phase 
transitions is still a subject of intense research. This is particularly true of homogeneous first-order phase 
transitions, i.e. phase transitions originating in the bulk without nucleation at surfaces or extrinsic defects. 
Recently, molecular dynamics calculations with ever increasing sophistication allowed detailed insights into the 
atomic mechanisms triggering homogeneous melting^"^. Experiments providing similarly detailed atomic-scale 
information about the very beginning of the nucleation process, however, are scarce^. Some of the calculations 
suggest that spatial separation of defect pairs is a crucial step in early nucleation^'^. Specifically, for the homo- 
geneous melting of ice a defect pair separation mechanism has been proposed^ which substantially reduces the 
free-energy barrier as compared to the purely statistical aggregation of defects until critical melting nuclei are 
formed. In the present study we investigate experimentally a homogeneous phase transition occurring in a two- 
dimensional system. Observation of the phase transition with atomic resolution reveals that it, too, proceeds via a 
defect separation mechanism. It can be shown that the process is analogous to phase transitions in spin systems. 
This finding suggests that separation of defect pairs is a common mechanism for reducing the nucleation barrier 
thus acting as a trigger for phase transitions in a wide variety of systems. Close to the phase transition temperature 
Tph defect pairs are formed, split, and the propagating moieties convert one phase into the other, much like sliders 
in a zipper mechanism. At temperatures different from Tph the defect pairs may still be formed, but their 
separation is hindered by a string interaction which thus suppresses the phase transition. Since the string 
interaction decreases as the dimensionality of the system is reduced, quasi- critical fluctuations can occur within 
a temperature range which increases with decreasing dimensionality. 

Results 

The phase transition. Here we investigate an order-order phase transition in the Br/Pt(110) adsorbate system at 
half-monolayer (ML) coverage from a (2 X 1) rectangular phase at low temperature to a c(2 X 2) checkerboard 
structure at room temperature. Previous analysis yielded the following results^"^: In both ordered phases every 
second short bridge site on the substrate close-packed rows is occupied, but with either zero phase difference from 
row to row ((2X1) structure) or with a phase difference of n (c(2 X 2); see Figure 1). The energetic preference for 
the (2X1) structure at low temperature is due to a coupling between charge density wave fluctuations associated 
with a periodic lattice distortion (CDW/PLD) and the Br adlayer. At low temperature the Bromine atoms lock 
into the CDW/PLD pattern and form a (2 X 1) structure. Since in first approximation the coupling is linear in the 
order parameter of the CDW, this has the effect of an external field stabilising the CDW/PLD. At room tem- 
perature, in contrast, electronic and lattice entropy lead to the disappearance ("melting") of the CDW/PLD and 
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Figure 1 | Ordered structures of 0.5 ML Bromine on Pt(llO): (a) Upper 
left and lower right: checkerboard structure with a c(2 X 2) unit cell 
containing 2 atoms; lower left and upper right: rectangular structure with a 
(2X1) unit cell containing 1 atom, (b) (3 X 2) unit cell containing 3 atoms 
embedded into c(2 X 2) and (2X1) domains. Black dotted rectangles: unit 
cells. Also shown are possible mechanisms of checkerboard ^ rectangular 
conversion, i.e. by travelling solitons (turquoise ellipse) or anti-solitons 
(red ellipse), and via a concerted place exchange involving a (3 X 2) unit 
cell. 

the repulsive energy in the adlayer is minimised by forming a per- 
fectly long-range ordered quasi-hexagonal c(2 X 2) configuration on 
a flat substrate^'^ 

Fluctuations. According to the Landau rules phase transitions 
between ordered states with different order parameters are always 
first order transitions. Nevertheless we observe in the present case 
quasi- critical fluctuations at 50 K indicating an anomalously low 
barrier between the two ordered states. The fluctuations are 
illustrated in Figure 2 showing an STM image of a Pt(llO) surface 
covered with 0.5 ML of Br. The preparation shown here contains 
defects which allow a unique identification of the areas observed. 
Clearly, an inter- conversion of domains occurs in both directions. 
Note that in Figure 2 the length scale is set by the average defect- 
defect distance, but on a defect- free surface there is no preferred 
domain size as expected for quasi- critical behaviour. 

The defect-pair mechanism. In order to investigate the mechanism 
for the fluctuations on an atomic scale, we first consider a 
preparation which is virtually free of extrinsic point defects, such 
as shown in Figure 3a. If one zooms in to examine the transition 




Figure 2 | STM images (12.6 X 12.6 nm') of a defective Br/Pt(110) 
surface at 50 K recorded in succession within a time interval of about 

40 s. The Br coverage is precisely 0.5 ML. White ellipses mark identical 
areas where checkerboard structures were converted into rectangular ones 
and vice versa. 

region between (2X1) and c(2 X 2) domains one usually observes a 
reduced and blurred contrast, sometimes even a loss of atomic 
resolution, although within the domains the atoms are well 
resolved. The blurring is attributable to the mobility of atoms at 
the domain boundaries. On careful examination, however, the 
transition from local (2X1) into c(2 X 2) domains and vice versa 
is revealed to be caused by either a compression (Figure 3a) or a 
dilution (Figure 3c) of Br atoms along the [lIO] direction (the 
close-packed row direction of the substrate). Figure 3b shows how 
such intrinsic defects are generated in an otherwise perfectly ordered 
domain by a density fluctuation along the [iTO] direction. In the row 
marked by white lines a displacement of a Br atom from its regular 
short bridge adsorption site towards one of its neighbours is found. 
This creates a dilution followed by a compression in the Br atom 
density along the row. The defect pair can be considered as a 
combination of soliton and anti-soliton, respectively, since each 
component introduces an equally large phase difference tt, but 
with opposite sign, between neighbouring rows. If soliton and anti- 
soliton are formed in a c(2 X 2) domain and separate from each 
other, then, like sliders in a zipper, they leave a stripe of (2 X 1) phase 
with a width of three adjacent rows in their wake. Conversely, if they 
are formed in a (2 X 1) domain and separate, they create a c(2 X 2) 
domain. For a soliton (dilution) or an anti-sohton (compression) to 
travel along the rows, only successive single-atom jumps are required 
as illustrated in Figure la. 

Extrinsic defects. The presence of extrinsic defects changes the 
mechanism for the phase transition substantially. As seen in 




Figure 3 | Checkerboard-rectangular domain boundaries in a 0.5 ML Br/Pt(l 10) preparation virtually free of extrinsic point-defects, (a) At the centre 
of the white ellipse a c(2 X 2) ^ (2 X 1) transition occurs via a compression of the Br-Br distances along the [1 10] direction, (b) In the row marked by 
white lines a Br atom is seen to be momentarily displaced towards its top neighbour creating a density fluctuation which consists of a compression 
(soliton) paired with a dilution (anti-soliton). (c) Transitions between c(2 X 2) and (2X1) structures via anti-solitons. Black circles mark Br atoms with 
the regular succession of two Pt nearest-neighbour distances. White circles mark Br atoms with slightly larger distances thus accumulating a phase shift of 
71 from bottom to top. 
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Figure 4 | Schematic of the zipper mechanism converting a checkerboard 
into a rectangular phase as observed in the present surface system. Red 

and green arrows, respectively, designate the pseudospin, which is defined 
as +1/2 for occupied and —1/2 for unoccupied short-bridge sites. White 
arrows indicate the hopping of Br atoms. White ellipses mark solitons and 
anti-solitons (double particles and double vacancies, respectively). 



Figure 2, transitions between the two domains with different order 
parameter are in this case always associated with local (3 X 2) unit 
meshes (see Figure lb) where in addition to short-bridge sites one 
long bridge site per unit cell is occupied. On average the (3 X 2) 
arrangement is only 18 meV per Br atom less favourable than the c(2 
X 2) order, but - as for the (2X1) phase - only if the substrate is 
allowed to relax into a particular buckling pattern. Local nucleation 
of the (3 X 2) pattern is therefore associated with an energetic barrier 
higher than the value above, since the latter refers to a global periodic 
elastic distortion of the substrate. The nucleation barrier is 
apparently significantly higher than for the (2X1) buckling. In 
the presence of extrinsic point defects, however, the barrier is 
reduced and spontaneous nucleation of (3 X 2) structural elements 
is observed. Once they are nucleated, the (3 X 2) structure units play 
a dominant role in the inter- conversion of the c(2 X 2) and the (2 X 
1) phase, as can be concluded from a detailed inspection of Figure 2. 
The mechanism is sketched in Figure lb: The (3X2) unit mesh can 
propagate along the [lIO] direction via a concerted hopping of three 
neighbouring atoms in a Ring-a-Ring-o' -Roses round dance fashion. 
As the (3 X 2) unit mesh is shifted along the close-packed row 
direction the inter- conversion of the c(2 X 2) and the (2X1) 
phase proceeds. 

Discussion 

Introducing the occupation numbers ±1/2 for occupied and unoc- 
cupied short-bridge sites, respectively, i.e. assigning a pseudospin to 
each short bridge site, allows mapping of the present problem onto a 
spin system. This is shown in Figure 4 for a c(2 X 2) ^ (2 X 1) 
conversion. A density fluctuation translates into a pair of flipped 
spins (Figure 4a). This creates a soliton anti-soliton defect pair 
(Figure 4b). The defect moieties separate and travel independently 
through the system thereby creating via a zipper mechanism a stripe 
of the complementary phase between them. Far away from the phase 
transition temperature Tph appearance of the "wrong" phase raises 
the free energy of the system. This amounts to an effective attractive 
interaction between the defects and is known as string interaction. As 
a result, soliton and anti-soliton tend to recombine which restores 
the original order. Close to Tp^ the free energy of the two phases is 
balanced, the string interaction disappears, and the soliton pair 



becomes unbound. The travelling solitons result in quasi- critical 
phase fluctuations. 

Figure 4 emphasizes the similarity of the present mechanism to 
phase conversions in strongly correlated systems. For instance, 
removal of an electron from an antiferromagnetic system creates a 
spin (spinon) and charge (holon) excitation which can also be con- 
sidered as a soliton anti-soliton defect pair. In a strictly one-dimen- 
sional (ID) system (Tomonaga-Luttinger-Liquid; TLL), there is no 
string interaction. Thus spinon and holon may always freely separate 
and travel independently along the chain^°. In a higher- dimensional 
antiferromagnet the string interaction prevents spin and charge 
separation, much as in the present case far from Tp^ the soliton 
anti-soliton separation is prevented. 

The analogy between spin and particle chains made so far is form- 
ally reflected in the equivalence between the Heisenberg and the 
Bose-Hubbard Hamiltonian^"'^^ respectively. Of course, the pro- 
pagation of the quasi-particles in the present case is not coherent 
in contrast to these models, since the Br atoms are coupled to the 
phonon bath of the substrate. On the other hand, the hopping of the 
Br atoms is not simply a thermal diffusion process determined by 
the diffusion barriers along the chain on a rigid substrate. In that case 
one should expect considerable dynamics at room temperature, 
while the measurements reveal an essentially static c(2 X 2) pattern 
at 300 K^. In contrast, the hopping is apparently associated with 
thermally driven CDW/PLD fluctuations occurring at 50 K. We have 
previously proposed an extended Landau model which qualitatively 
describes the energetics of the phase transition^: 

F= ^/s/S/+i+ y^^m(5/ + 5/+i)+ ^^m^+ ^m^ (1) 

Here, F is a free- energy density, m is the order parameter of the 
CDW/PLD in the substrate and the last two terms are the customary 
Landau expansion of the free energy in terms of the order parameter 
for the CDW phase transition. / > 0 is the repulsive energy between 
Br atoms on neighbouring short-bridge sites across the rows (i.e. nnn 
short-bridge sites). S/, S/+i are the corresponding occupation num- 
bers (pseudo spins), as described above. ^ < 0 is a coupling constant 
accounting for the coupling between the adsorbate atoms and the 
CDW in the substrate. Due to the adsorbate-substrate coupling^ the 
Br atoms are dragged along^ with the CDW fluctuations, i.e. fluctua- 
tions of the order parameter m which occur around the CDW phase 
transition temperature. With increasing temperature the fluctua- 
tions become faster and the atomic resolution is lost^ As T exceeds 
200 K, however, the CDW/PLD fluctuations die out, since m ^ 0 
and, owing to the repulsive Br-Br interaction represented by the first 
term in equation (1), a static c(2 X 2) Br pattern develops. 
Association of the Br mobility with CDW fluctuations is supported 
by DFT calculations, showing that on a rigid substrate there exists a 
substantial barrier of the order of 1.2 eV for the occupation of neigh- 
bouring short-bridge sites^^'^^. 

In summary, an atomically resolved observation of a first- order 
phase transition in a 2D system reveals a defect-pair formation as the 
crucial element initiating the transformation. Close to the transition 
temperature the composite defect can decay into a compression 
(soliton) and a dilution (anti-soliton) travelling independently 
through the 2D crystal. In a kind of zipper mechanism they create 
a stripe of the altered phase in their wake. Far from the transition 
temperature the defect pair is bound by an effective string inter- 
action. Since the string interaction decreases with dimensionality, a 
strongly anisotropic, quasi- ID system such as the present one can 
exhibit quasi- critical fluctuations in a comparatively wide temper- 
ature range around Tph even though the phase transition described 
by equation (1) is weakly first-order. The soliton anti-soliton zipper 
mechanism partially circumvents the conventional nucleation bar- 
rier. The mechanism reported here parallels similar processes in 
correlated systems and also resembles a recently postulated mech- 
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anism for ice melting^. Therefore it is apparently of general relevance 
and one may expect similar zipper mechanisms to be uncovered in a 
wide variety of systems, once phase transitions can be observed with 
atomic resolution. 

Methods 

Experimental. Our model system has been described previously^'^^'^^. Almost defect- 
free c(2 X 2)-Br/Pt(110) surface preparations were obtained in a preparation 
chamber with base pressure of 7 X 10"^^ mbar. Deposition of a few monolayers of Br 
onto the clean ( 1 X 2)-mr reconstructed Pt( 110) surface was followed by annealing to 
about 780 K. This results in partial Br desorption, lifting of the Pt(llO) missing- row 
reconstruction and evolution of a long-range ordered c(2 X 2) structure as verified by 
low-energy electron diffraction (LEED). For STM investigations the sample was 
freshly prepared (no LEED check, to avoid the introduction of defects) and 
transferred into the STM chamber with a base pressure of ~3 X 10"^^ mbar. 
Topographic images were recorded with an Omicron VT-STM featuring a He flow- 
cryostat. Defective surfaces were prepared by using a marginally lower annealing 
temperature after Br deposition. One class of defects observed after such a 
preparation consists of single dark objects which we attribute to Pt atom defects 
remaining after the lifting of the ( 1 X 2)-mr reconstruction. Other defects appeared as 
bright objects in the STM image. Some of them seemed to be related to the H2O 
background pressure (total background pressure 10"^^ mbar range) and showed a 
tendency for island formation. We tentatively attributed them to H2O adsorbing after 
prolonged cooling. 

Calculations. DFT calculations have been performed with the Vienna Ab-initio 
Simulation Package (VASP)^'*'^^ using the projector augmented wave method^^ for 
periodically repeated slabs of varying substrate thickness and adsorption geometries 
separated by a vacuum region of minimal 14 A. In order to assess the effects of 
substrate thickness, adsorption geometry, and the particular functional used, we 
performed all calculations for three different (semi) -local functional, LDA^^, PBE^^, 
and PBEsoP'', which yield different Pt substrate lattice constants a, with aLOA (3.91 A) 
slightly smaller, apBE (3.97A) larger and apBEsoi reproducing the experimental value of 
3. 92 A. The influence of the finite substrate thickness was judged by performing all 
calculations on a 9 layer thick Pt slab with Br adsorbed on one surface, and on 1 1 layer 
and 15 layer thick Pt substrates with Br adsorbed on both surfaces. For the latter two 
configurations the interactions between the two Br layers across the substrate was 
studied by using both a mirror-plane and inversion geometry configuration of the two 
Br layers. All atomic positions were fully relaxed until the largest force was below 
0.01 eV/A. Averaging the energy difference between the c(2 X 2), and (3 X 2) 
geometries over all functional and slab models leads to a mean square deviation of 
17 meV. This spread is rather caused by the different functional than by the finite 
slab models, for which it amounts to only 7 meV. 
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